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Abstract. We study the changes in physical and chemical conditions during the early stages of collapse of a 
pre-protostellar core, starting from initial conditions appropriate to a dense molecular cloud and proceeding 
to the “completely depleted” limit. We allow for molecular desorption from the grain surfaces and follow the 
evolution of the ionization degree and the ionic composition as functions of time and density. The timescale for 
collapse is treated as a parameter and taken equal to either the free-fall or the ambipolar diffusion time. The 
processes of freeze-out on to the dust grains and of coagulation of the grains were treated simultaneously with 
the chemical evolution of the medium in the course of its collapse. When proceeding at close to its maximum rate, 
coagulation has important consequences for the degree of ionization and the ionic composition of the medium, but 
its effect on the freeze-out of the neutral species is modest. An innovation of our study is to calculate the grain 
charge distribution; this is done in parallel with the chemistry and the dynamics. The grain charge distribution 
is significant because ions recombine predominantly on the surfaces of negatively charged grains. We have 
also attempted to reproduce with our models the observational result that nitrogen-containing species, such as 
NHg and N 2 H+, remain in the gas phase at densities for which CO and other C-containing molecules appear to 
have frozen on to grain surfaces. We conclude that recent measurements of the adsorption energies of N 2 and CO 
invalidate the interpretation of these observations in terms of the relative volatilities of N 2 and CO. We consider an 
alternative explanation, in terms of low sticking coefficients for either molecular or atomic N; but this hypothesis 
requires experimental confirmation. We find that, irrespective of the nitrogen chemistry, the main gas phase ion 
is either or (and its deuterated isotopes) at densities above 10® cm“®; whether H'*' or predominates 
depends sensitively on the rate of increase in grain size (decrease in grain surface area per unit volume of gas) 
during core contraction. Our calculations show that will predominate if grain coagulation proceeds at close 
to its maximum rate, and otherwise. 
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1. Introduction 

The contraction of an interstellar “core” to form a proto¬ 
star is a process of great current interest. It is clear that 
this contraction is accompanied by considerable changes 
both in the characteristics of the grains associated with 
the core and the chemical composition of the gas. In par¬ 
ticular, most molecules containing heavy elements con¬ 
dense on to grain surfaces (Bacmann et al. 2002, Tafalla 
et al. 2004). For CO and several other species, this process 
occurs at densities uh ^ 310^ cm“^. Moreover, there is 
evidence in several objects for a “molecular hole” at den¬ 
sities above 10® cm“®, where material containing C,N,0 
is in solid form and only molecules and ions composed of 
H and D remain in the gas phase. 
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In two earlier articles (Walmsley et al. 2004, Flower et 
al. 2004), we have considered the characteristics of such 
“completely depleted” regions and have concluded that 
the degree of ionization of the gas and the relative abun¬ 
dances of the major ions (H+, Hg , H 2 D+, D 2 H+, Dg ) are 
sensitive to the size of the negatively charged grains, with 
which H+ ions recombine. Recent observations of H 2 D“'' 
and D 2 H+ (Caselli et al. 2003, Vastel et al. 2004) are 
compatible with the existence of a zone lacking heavy ele¬ 
ments. However, detection of these ions (H 2 D+ and D 2 H+) 
seems to be rare, suggesting that objects containing such 
a depleted nucleus are in the ultimate phase of evolution, 
prior to the formation of a protostar. It seems relevant also 
that the observations suggest that there is an intermedi¬ 
ate density range (10® ^ uq ^ 10® cm“®) where molecules 
containing nitrogen, such as NHg and N 2 H+, remain in the 
gas phase when carbon- and oxygen-containing species 
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are already depleted out. This effect of differential deple¬ 
tion was believed to be due to the relatively high volatility 
of N 2 (see Bergin & Langer 1997), which remains in the gas 
phase when CO, for example, has condensed out. However, 
recent measurements of the adsorption energies of CO and 
N 2 , which are found to differ by less than 10% (Oberg et 
al. 2004), do not support this interpretation. Accordingly, 
we attempt to account for the differential depletion by 
varying other (uncertain) parameters of the model. 

An important issue is the extent to which grain co¬ 
agulation takes place during contraction, thereby influ¬ 
encing the timescale for molecular depletion. Studies of 
this process (e.g. Ossenkopf 1993, Ossenkopf & Henning 
1994, Suttner & Yorke 2001) suggest that, at densities 
10® ^ ^ 10^ cm“®, grains grow in size, although not 

to the same extent as in protoplanetary disks, where the 
density is higher. The grain opacity per H nucleus, de¬ 
fined here as (ng/nH)7rag, where rig is the number density 
of grains of radius Og and nn = n(H)-|- 2 n(H 2 ), is expected 
to decrease but not sufficiently to prevent depletion on to 
grain surfaces (whose rate is proportional to the same pa¬ 
rameter); this is consistent with what is known about the 
mm-sub-mm opacity of prestellar cores (e.g. Bianchi et 
al. 2003, Kramer et al. 2003). 

With the above points in mind, we have undertaken 
calculations of the evolution of the chemical composition 
of such a contracting core. Our work is similar in some re¬ 
spects, but differs in others, to studies of the surface and 
gas-phase chemistries in a contracting core by Aikawa et 
al. (2001, 2003, 2005), as well as the time-dependent cal¬ 
culations of Roberts et al. (2003) and the recent study 
of Bonnor-Ebert spheres by Lee et al. (2004). Our model 
incorporates a less detailed treatment of the gas-phase 
chemistry, and we do not attempt to treat surface reac¬ 
tions explicitly. Instead, we have emphasized the evolu¬ 
tion of grain charge and size with time and the conse¬ 
quences for the ion chemistry, which we perceive to be a 
key aspect of the problem. The grain size (through the 
grain surface area) influences the relative abundances of 
H+ and H^^, which are likely to be the major ions during 
the later phases of core contraction. Our models bring out 
the important role played under some circumstances by 
H+, whose significance appears to have been neglected - 
or, at least, not emphasized - in earlier studies. 

In Section 2 of this article, we summarize our formula¬ 
tion of the problem. In Section 3, we consider the critical 
input parameters of the model. Our main results are given 
in Section 4 and our concluding remarks in Section 5. 

2. Formulation 

We consider the evolution of a contracting sphere (“core”) 
of gas and dust and, for simplicity, we use equations appli¬ 
cable to free-fall, homologous, isothermal collapse. This 
procedure allows us to parametrize the dynamical be¬ 
haviour in terms of a collapse timescale, Tc, which we vary 
between a minimum value, taken to be the free-fall time, 
rg, and a maximum value, defined by the ambipolar dif¬ 


fusion timescale. Tad (see below), both evaluated accord¬ 
ing to the initial conditions in the medium. In practice, 
the ambipolar timescale is about an order of magnitude 
larger than the free-fall time. Thus, Tc spans timescales 
between the magnetically subcritical case, in which the 
magnetic flux suffices to prevent collapse, and the mag¬ 
netically supercritical case, where collapse can take place 
immediately. 

The equation which determines the variation of the 
radius R of the core with time t is 

1 di? _ 1 dec ^ ^ ^ ('ll 

R dt X dt 2 tcX / 


where x = R/Rq < \ and Rq is the initial radius of 
the condensation. The timescale for free-fall collapse is 
(Spitzer 1978) 

1 . 

2 
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where po is the initial mass per unit volume. In spherical 
collapse, the mass density p varies with time t according 
to 


(3) 

p dt X dt 

In practice, the collapse is described by equation CJ, with 
the timescale Tc taken to be either rg or the ambipolar 
diffusion timescale. Tad, to be defined below. As already 
noted. Tad is about an order of magnitude larger than Tg. 

The medium contains neutral gas, positive ions and 
electrons; it also contains neutral and (mainly negatively) 
charged grains. The number density, n, of the gas-phase 
neutrals is determined by 

1 dn 3 dx ^ A/" 

n dt X dt n 


where TV is the rate per unit volume at which neutrals are 
created in the gas phase by chemical reactions. Similarly, 
for the positive ions in the gas phase 

1 dn+ 3 dx TV+ . , 

-^ = —+ — 

n+ dr X dr n+ 


Analogous equations are solved for each of the chemical 
species which intervene in the chemistry. 

Central to our approach, and distinguishing it from 
previous work, is the self-consistent calculation of the 
change in the grain cross section which occurs owing to 
coagulation and to freeze-out from the gas phase. The 
grains are assumed to have a unique radius, determined 
initially by the graimgas mass ratio and the mass den¬ 
sity of the grain material. Subsequently, the grain radius 
varies owing to the freeze-out of gas-phase constituents 
on to the surfaces of the grains and to grain coagulation. 
The rate of variation of the grain number density is given 

by 

1 drig 3 dx ^ TVg 

rig dt X dt rig 
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where rig is the grain number density. The first term on 
the right hand side of this equation represents the rate 
of increase of the number density owing to the compres¬ 
sion associated with the (spherical) collapse; the second 
represents the rate of decrease of the grain number den¬ 
sity owing to coagulation. When two grains stick together, 
only one grain remains, and hence the rate per unit vol¬ 
ume at which the grains are removed through coagulation 
is A4 = — {av)n‘^. The method used to determine the rate 
constant (av) is presented in Appendix A. We have as¬ 
sumed the grains to be spherical, which minimizes their 
surface area and effective cross section for a given volume. 
Clearly, non-spherical grains would have a larger cross 
section. 

In addition to the free-fall timescale [equation lO] , we 
evaluate the ambipolar diffusion timescale from 


'^ad — 


^ rii 


mimn 


7rGm^ - -r' rin rrii + 


(7) 


(cf. Flower & Pineau des Forets 2003a) for both positive 
ions and charged grains (in which case rrii » m-n) and 
adopt the larger of the two values. We use the expression 
of Osterbrock (1961) for the rate coefficient for momen¬ 
tum transfer in a collision between an ion and a neutral 
particle, 

((Tu)in = 2.4l7r(e^a„/min)^ (8) 


where niin = WiTTin/ imi + m-o) is the reduced mass and On 
is the polarizability of the neutral. For collisions between 
neutrals and charged grains, we take a = TrOg and v = 
{SkBTn/TTrrin) 2 , where Tn is the temperature of the neutral 
gas. 


3. Model 

Adopting the initial conditions in steady-state, the con¬ 
densation was allowed to undergo isothermal collapse, as 
described in Section 2. The dynamical and chemical rate 
equations were solved simultaneously. The grain number 
density was determined by the competing effects of com¬ 
pression, arising from the collapse, and coagulation, which 
tends to reduce the total number of grains [cf. equation 0 
above]. The grain radius was calculated self-consistently 
as a function of time, t, from the grain number density and 
the total mass of material currently in the solid phase, al¬ 
lowing for the continuing adsorption of heavy species on to 
the grains. In practice, adsorption has only a minor effect 
on Og, as a mantle of ices is assumed to exist before the 
onset of collapse; coagulation is a much more significant 
mechanism for increasing Og. 

3.1. Initial conditions 

The initial distribution of the elemental abundances across 
the gas phase and the solid phase (i.e. grains) is given 
in Table ^ The elemental abundances derive from the 


Table 1. The elemental abundances and the initial repar¬ 
tition between the gas phase and the solid phase: based 
on the studies of Anders & Grevesse (1989), Savage & 
Sembach (1996), Gibb et al. (2000), and Sofia & Meyer 
(2001). ‘Depletion’ denotes the fraction of the element in 
the solid phase (grains). Numbers in parentheses are pow¬ 
ers of 10. 


Element 

Fractional abundance 

Gas phase 

Solid phase 

Depletion 

H 

1.00 

1.00 

3.95(-4) 

0.0 

He 

l.OO(-l) 

l.OO(-l) 


0.0 

C 

3.55H) 

8.27(-5) 

2.72(-4) 

0.77 

N 

7.94(-5) 

6.39(-5) 

1.55(-5) 

0.20 

0 

4.42(-4) 

1.24(-4) 

3.18^4) 

0.72 

Mg 

3.70(-5) 


3.70^5) 

1.0 

Si 

3.37^5) 


3.37^5) 

1.0 

S 

1.86^5) 

0.60(-6) 

1.80^5) 

0.97 

Fe 

3.23^5) 

1.50(-9) 

3.23(-5) 

1.0 


study of Anders & Grevesse (1989), and the depletions 
from Savage & Sembach (1996), Gibb et al. (2000), and 
Sofia & Meyer (2001), with the exception of sulphur, for 
which the depletion depletion from the gas phase within 
molecular clouds is difficult to establish observationally. 
We have assumed that sulphur is depleted from the gas 
phase by a factor of 30 (cf. Wakelam et al. 2004). Given 
the distribution in Tabled the initial value of the mass of 
dust to the mass of gas is 0.0094. 

The computer code was run to determine the steady- 
state abundances of the chemical species in the gas phase 
(neglecting any further depletion) by the simple expedient 
of specifying a timescale for collapse, Tc, in equation o 
which is much greater than the chemical equilibrium 
timescale; 92 gas-phase species and almost 800 chemical 
reactions were included in the model. The gas-phase 
chemistry included ion-neutral, neutral-neutral, and 
dissociative recombination reactions involving species 
containing H, He, C, N, O and S. The branching ratio 
for the dissociative recombination of N 2 H+, to N or N 2 , 
was taken from the recent measurements of Geppert et 
al. (2004). Ionization by cosmic rays, both direct and in¬ 
direct (by the electrons released), was taken into account. 
Adsorption of neutral species on to grains, and desorp¬ 
tion of the hydrogenated products of the grain-surface 
chemistry following cosmic ray heating of the grains, 
were also incorporated (see below). The rate of cosmic 
ray ionization of hydrogen was taken to be ^ = 1 x 10“^^ 
s“^. The list of species and the complete set of chem¬ 
ical reactions included in the model are available from 
http://massey.dur.ac.uk/drf/protostellar/speciesjchemistry 

The initial value of the grain radius was taken to be 
Og = 0.05 ^m, which yields a grain opacity, defined here 
as rigTrOg, of 1.6 x lO^^^nn cm“^, assuming the density of 
the grain material to be 2 g cm“^. Adopting the grain- 
size distribution of Mathis et al. (1977), with limits of 
0.01 < Og < 0.3 fim to the distribution, and assuming the 
canonical value of 0.01 for the dust:gas mass ratio, one 
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obtains the same value of Ugira^. [We note that the value 
of Og required to reproduce a given value of rigTrOg depends 
only on the adopted size distribution and its limits, and 
not on the density of the grain material nor the dust:gas 
mass ratio: see Appendix B.] 


3.2. Grain charge distribution 

In view of the significance of grains for the state of charge 
of the medium, the grain charge distribution was com¬ 
puted in parallel with the chemical composition, allowing 
for electron attachment and detachment and electron-ion 
recombination on grain surfaces (cf. Flower & Pineau des 
Forets 2003b, Section 2.2). Singly charged and neutral 
grains were considered. In the case of collisions betwen 
charged particles and oppositely charged grains, the ‘fo¬ 
cusing’ effect of the Coulomb attraction was included by 
means of equation (3.4) of Draine & Sutin (1987). In the 
case of collisions between charged particles and neutral 
grains, the attraction arising from the polarization of the 
grain by the incoming particle was included by means of 
their equation (3.3). For example, the rate per unit volume 
of electron attachment to neutral grains was evaluated as 
nengoTTOg)! -I- {ire'^/2agkBTe)^'^][8kBTe/ (7rme)]°'®Se, where 
rie is the electron density, rigo is the density of the neutral 
grains, Tg is the electron temperature (taken equal to the 
gas temperature, T = 10 K), nie is the electron mass, and 
Se = 0.5 is the electron sticking coefficient. The initial 
value of the total grain number density, rig, was calcu¬ 
lated self-consistently from the elemental depletions, the 
specified grain radius, and the mean density of the grain 
material, taken to be 2 g cm“^ at all times. The grain 
charge distribution is not sensitive to the assumed grain 
size, as the dominant charge-changing processes - elec¬ 
tron attachment to neutral grains and ion recombination 
with negatively charged grains - scale approximately with 

3 

a| and with Og, respectively, under the conditions which 
are relevant here. 


3.3. Adsorption and desorption 

The rate coefficient for adsorption of neutral species on to 
the grains was taken to be 7 rag[ 8 A:BT/( 7 rmn)]°'®S'n, where 
rrin is the mass of the neutral. The ‘sticking’ coefficient 
= 1, unless stated otherwise; we assumed that He re¬ 
mained in the gas phase and that H 2 was returned to the 
gas phase after forming on the grains. The geometrical 
cross section of the grain is adopted in this expression for 
the rate coefficient. Although the potential arising from 
the polarization of the grain by a charged particle can 
be significant (see Section 3.2), the potential due to the 
polarization of a neutral by a charged grain is negligible; 
this is because the polarizability of a grain is orders of 
magnitude larger than the polarizabilities of the neutral 
gas-phase species, owing to the much larger size of the 
grain. 


The recombination of positive ions with electrons on 
the surfaces of negatively charged grains releases the elec¬ 
tron binding energy, which is typically a few eV and 
large compared with the adsorption energies of the neutral 
species. Accordingly, we have assumed that the neutrals 
which are formed in this way are returned immediately to 
the gas phase. The attachment of ions to neutral grains 
is a relatively minor process, owing to the predominance 
of negatively charged grains (see Section 4.3 below) and 
to Coulomb focusing in collisions of ions with negatively 
charged grains. Consequently, we have neglected the ad¬ 
sorption of ions on to grains when computing the rate of 
freeze-out of gas-phase species. 

With the exception of H 2 , which desorbs following its 
formation, desorption of molecules was assumed to be con¬ 
trolled by the heating of grains which occurs owing to cos¬ 
mic ray impact (Leger et al. 1985). The rate of desorption 
of species i per unit volume per unit time from the grains 
is 


Tt- 

X^UgA-Kal 7 exp[ 
Z^i ’h 


_^^ads _ ^ads^ 

^max 

g 


where is the fractional abundance of species i 

on grains, is the adsorption energy of the species 

(Aikawa et al. 1996, 1997, 2001), which relate to mea¬ 
surements on pure ices; = 70 K is the maximum 

temperature reached by the grains following cosmic ray 
impact (Hasegawa & Herbst 1993). We adopted 7 = 70 
cm“^ s“^, which yields a rate of desorption for CO which 
is consistent with the value for ‘spot’ heating, as derived 
by Leger et al. (1985); the heating arises from the impact 
of heavy cosmic rays, particularly Fe nuclei. In the case of 
small grains, ‘whole-grain’ heating may be more impor¬ 
tant than ‘spot’ heating, and the corresponding value of 
7 would be larger than adopted here. 


3.4. Coagulation 

Grain coagulation was treated as described in Appendix 
A. The important parameters are the local turbulent speed 
of the gas, which is assumed to determine the relative 
speed of the grains, and the critical speed Wcrit(ag) below 
which grains are assumed to coagulate on every collision. 
Following Tafalla et al. (2004), who observed L1498 and 
L1517B (see also Barranco & Goodman 1998), we adopted 
the turbulent line-width as an upper limit to the FWHM 
of the Gaussian velocity distribution, i.e. Av ^ 0.12 km 
s“^, which is approximately Cs/2 at T = 10 K, where Cg is 
the adiabatic sound speed. This value may include contri¬ 
butions from motions on scales larger than the mean-free 
path for grain-grain collisions and hence places an up¬ 
per limit on Av. Accordingly, we have taken Av = 
as our reference value, but we have also run models with 
Av = Cs/2 and Av = Cg/lO. 

The critical speed for coagulation, Vcrit) was taken from 
Chokshi et al. (1993), unless otherwise stated. We adopted 
their results for the case of ‘ice’; this seems the most ap- 



D.R. Flower et al.: Freeze-out and coagulation in pre-protostellar collapse 


5 


propriate of the three cases (‘ice’, ‘graphite’, ‘quartz’) con¬ 
sidered by Chokshi et ah, as the grains already have ice 
mantles at the start of the calculation. Accordingly, 

fcrit(ag) = 0.1 Og®/® (9) 

when Ucrit is in cm s“^ and Cg is in cm. We have assumed 
that the grains remain spherical throughout. This assump¬ 
tion minimizes the grain radius, and hence the grain opac¬ 
ity per H nucleus, (ng/nH)7rag, as the collapse proceeds. 
We assume also that the density of the grain material 
is unchanged by coagulation. The critical speed given by 
equation is greater than the estimate of Suttner & 
Yorke (2001) but approximately 4 times less than the mea¬ 
surement of Poppe & Blum (1997). We have run some 
comparison calculations in which Ucrit was taken to be 4 
times greater than is given by equation ©■ 

4. Results 

We consider the case of a condensation of gas and dust 
with an initial number density tt-h = n.(H) -|- 2 n(H 2 ) = 10^ 
cm“® and kinetic temperature T = 10 K. We recall that 
the rate of cosmic ray ionization of hydrogen was taken to 
be C = 1 X 10“^^ s“^. 

4.1. Freeze-out 

Fig. n illustrates the freeze-out of selected neutral species, 
including CO and N 2 , on to the grains and the evolution 
of the fractional abundances of the major ions, includ¬ 
ing H+, H;^ and HCO+. Results are shown for a refer¬ 
ence model in which coagulation was neglected and the 
timescale for collapse is determined by the free-fall time, 
Tc = Tff. The corresponding density profile, is also 

shown in Fig. ^ As will be seen below, collapse gets un¬ 
der way after a few times 10® yr, and the neutral species 
then freeze out rapidly on to the grains. Initially, HCO'*' 
is the most abundant molecular ion, but it is overtaken 
by and H+ as freeze-out proceeds. Ultimately, the de¬ 
gree of ionization of the gas is determined by under 
the conditions of this model. 

The adsorption energy of N 2 which was used to obtain 
the results in Fig. ^ was 790 K, which is to be compared 
with = 855 K (Oberg et al. 2004). A much larger 
difference between the adsorption energies of N 2 and CO 
is required for there to be a significant difference between 
the density-dependence of the fractional abundances of 
N 2 and CO, and hence of N 2 H+ and HCO'*'. Another pos¬ 
sible explanation of this differential effect will be consid¬ 
ered below, in the context of a discussion of the nitrogen 
chemistry. 

The effect of allowing the ambipolar diffusion timescale 
to regulate the collapse is shown in Fig.[2t- As we shall see 
below, the timescale for ambipolar diffusion is about an 
order of magnitude larger than the free-fall time, and so 
freeze-out occurs at lower values of the gas density, before 
the collapse is under way. Indeed, in this case, freeze-out 


occurs at densities which are lower than the values of nn ^ 
3 X 10^ cm“® which have been deduced from observations 
of CO (cf. Bacmann et al. 2002, Tafalla et al. 2004). Thus, 
the observational data pertaining to CO depletion favour a 
collapse timescale of order of the free-fall time. The ‘knee’ 
which is apparent in the curves arises at the fractional 
abundance at which cosmic ray desorption begins to be 
significant; this occurs marginally sooner for N 2 than for 
CO, owing to the smaller (by 65 K) adsorption energy of 
the former species. 

Including coagulation, with the expression of Chokshi 
et al. (1993) for Ucrit(ag) [equation 0], has very little ef¬ 
fect on the results, as may be seen by comparing Figs.Et 
and Eb: the ‘free-fall’ results in Fig. Eb (the upper pair 
of curves) are almost coincident with the full curves in 
Figs. Eb- Even when the (4 times larger) values of Ucrit(ag) 
are adopted, following Poppe & Blum (1997), which pos¬ 
sibly yields an upper limit to the rate of coagulation, its 
effect on the evolution of the fractional abundances of the 
neutral species is not large, as may be seen from Fig. Eb- 

4.2. Timescales 

In Fig. El we plot the grain radius, Ug, as a function of the 
gas density, nn • Results are shown for the limiting case of 
no coagulation (ucrit = 0) and for Ucrit(ag) = 0.1 flg 
[equation 0; Chokshi et al. (1993)] and Ucrit(ag) = 
0.4 (Poppe & Blum 1997); Av = Cs/4, where 

Cs = 0.243 km s“^ is the adiabatic sound speed. When 
equation 0 is adopted, the initial value of the critical 
speed for coagulation is Ucrit = 0.0262 km s“^, almost an 
order of magnitude smaller than Cg. 

Fig. El shows that, when no coagulation occurs, the 
grain radius, Og, is practically constant, as the contri¬ 
bution of continuing freeze-out to the size of the grain 
mantles is small; we have assumed that the grains have a 
mantle of ices before collapse begins. The rate of coagu¬ 
lation, and hence the rate of increase of the grain radius, 
Og, is seen to be sensitive to the critical speed, Ucrit(ag): 
increasing the constant of proportionality in equation 0 
from 0.1 to 0.4, which is in line with the measurements of 
Poppe & Blum (1997), has the effect of enhancing signif¬ 
icantly the degree of coagulation over the free-fall time. 
Reducing the turbulent width, Av, which is the FWHM of 
the Gaussian distribution [equation (IA.3ll ]. has a similar 
effect to increasing Ucrit- In our model, the grain opacity 
per H nucleus, (ng/nH)7rag, varies inversely with the grain 
radius, Og, and so the effect of coagulation is to reduce the 
grain opacity and hence increase the timescale for deple¬ 
tion. The assumption which we have made, that the grains 
remain spherical, minimizes their surface area (and hence 
their geometrical cross section) for a given density of the 
grain material, which we take to be constant during coag¬ 
ulation. Although these assumptions are not realistic (cf. 
Ossenkopf and Henning 1994), they allow a lower limit to 
be placed on the influence of coagulation on freeze-out. 
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The gravitational free-fall time, rg = 0.43 x 10® yr for 
our initial density njj = 10 ^ cm“®, is much smaller than 
the initial ambipolar diffusion timescales for both positive 
ions and charged grains, t‘™ = 3.4 x 10® yr and = 

3.6 X 10® yr. The timescale for freeze-out of the neutral 
species on to the grains - the inverse of the adsorption 
rate - is given by 

Tfo = [ngTrOgUthS'] ^ (10) 

where 


fth 


Trmn ) 


is the thermal speed, is the kinetic temperature and 
TOn is the mass of the species; S is the sticking coefficient. 
Taking Tn = 10 K and = 28 mn (corresponding to CO, 
as an example) and S' = 1 , we obtain Tfo = 2.2 x 10 ® yr 
initially, when Og = 0.05 /im and nn = 10^ cm“®. Thus, 
the timescale characterizing freeze-out is comparable with 
the free-fall timescale and much smaller than that for am¬ 
bipolar diffusion. The timescale for coagulation is 

Teg = [ngdTTOgUrel] ^ (11) 

where Urei is the relative collision speed. The maximum 
possible value of Urei is Ucrit [cf. equation @]; adopting 
Urei = Ucrit, we See that 


'^io '^crit 
- ^ ^ 


( 12 ) 


neutral species with ionization potentials less than that 
of H) decrease. This phenomenon is illustrated in Fig. 01 
where the fractional abundances of H+ and a selection of 
other ions are plotted as functions of the gas density, nn. 
The signihcance of coagulation may be seen by comparing 
Fig. 2Ji, in which coagulation is neglected, with Figs. QJ) 
and c, in which coagulation is included. 

The recombination of H'*' with free electrons is a ra¬ 
diative process and hence slow. Consequently, H+ ions 
recombine with electrons predominantly on the surfaces 
of negatively charged grains. On the other hand, the re¬ 
combination of with free electrons is a non-radiative 
(dissociative) process and intrinsically fast. Using the pa¬ 
rameters adopted in our model, the ratio R of the rate 
(s“^) of dissociative recombination of in the gas phase 
and its rate of recombination on the surfaces of the (pre¬ 
dominantly) negatively charged grains is 

i? 2 X lO^^XeUg 

where = nelriYi and Og is expressed in /rm. Thus, when 
ttg = 0.05 /im, R = 1 is attained when Xe « 2 x 10“®; then, 
Hg recombines at equal rates on the surfaces of negatively 
charged dust grains and in the gas phase. Referring to 
Fig. EK, we see that this point is reached when uh ~ 10 ® 
cm“®. At lower densities, dissociative recombination pre¬ 
dominates, whereas, at higher densities, recombination on 
the surfaces of negative grains predominates, and the frac¬ 
tional abundances of and H®" vary similarly with tih . 
The ratio of the fractional abundances of H;^ and H+ 

n(H+)/n(H+) « 18/(1 + i?) 


With the values of the parameters specihed above, we ob¬ 
tain Teg = 1.8 X 10® yr initially, which is approximately 
equal to Tfo. However, Teg tends to increase with time, ow¬ 
ing to the increase in Og and consequent decrease in Uerit- 
As we have assumed the maximum possible value of Urei 
when evaluating Teg, we conclude that freeze-out tends to 
occur more rapidly than coagulation, providing the stick¬ 
ing coefficient 5” « 1 . 

4.3. Fractional ionization 

The signihcance of ambipolar diffusion and, in particular, 
the role of elastic scattering between neutrals and charged 
grains in protostellar collapse have been recognized for 
some considerable time (Nakano & Umebayashi 1980; Shu 
et al. 1987; Ciolek & Basu 2000). Whether the timescale 
for collapse, Tc, is determined by ambipolar diffusion or by 
gravitational free-fall depends on the fractional ionization 
of the gas [cf. equation I?))] and the degree of charge of the 
grains. 

At early times, is the most abundant atomic ion, 
and its abundance is maintained by charge transfer of H+ 
to S; the source of H+ is cosmic ray ionization of hydro¬ 
gen. When freeze-out of the neutrals on to the grains 
occurs, the fractional abundance of H+ tends to rise as 
losses through charge transfer reactions (with S and other 


for the parameters adopted in our model; this expression 
agrees well with the numerical results presented in Fig.EJi 
in the high density limit, where i? —> 0. [We recall that, 
owing to the absence of deuterated isotopes in the present 
model, the number density of “Hg ” should be understood 
as n(H 3 ) -I- n(H 2 D+) -|- n(HDj) -|- 71 ( 03 ).] 

In the case where there is no coagulation, as in Fig. Hi, 
the grain opacity per H nucleus, (ng/nH) 7 rag, remains 
practically constant as the collapse proceeds; cf. Fig. Eb. 
The H®" ions are removed (by recombination with elec¬ 
trons on the surfaces of negatively charged grains) on a 
timescale which is inversely proportional to rigirag and 
hence to tih. For tih ^ 10® cm“®, ions are removed pre¬ 
dominantly by dissociative recombination with free elec¬ 
trons, on a timescale which is inversely proportional to 
approximately n^®. As a consequence, the fractional abun¬ 
dance of ions initially increases relative to H+ as tih 
increases, and the case of ‘no coagulation’ is characterized 
by n(H+) > n(H+) at high gas densities. On the other 
hand, when coagulation occurs (assuming Av = Cs/4 in 
Fig.0, {ng/nii)Trag decreases as nn increases (cf. Fig.Eb), 
and the rate of recombination of the H+ ions on grains is 
reduced. Owing to the higher fractional abundance of the 
free electrons, 7 i(H;]") < 77 (H+) at high gas densities when 
i'crit(ag) = 0.4 Og^^^. These results are summarized in 
Figs. ^-c. 
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The results in Fig.^Jshow that the ratio n(H 3 +)/n(H+) 
falls substantially when coagulation is introduced: grain 
growth favours H+. On the other hand, the fact that H 2 D+ 
and D 2 H''' have been observed in pre-protostellar objects 
(Caselli et al. 2003, Vastel et al. 2004) suggests that the 
deuterated forms of may be the dominant gas-phase 
ions. If so, the present calculations imply that coagula¬ 
tion is not significant in reducing the grain opacity per H 
nucleus during the early phases of collapse. 

The final panel (d) in Fig. 01 shows the effect of re¬ 
ducing the value of the sticking coefficient for molecular 
nitrogen to S'(N 2 ) = 0.1 and should be compared with 
those in panel (c), where 5 '(N 2 ) = 1 , the default value 
for all species which freeze out; the results in Fig. 0}! are 
considered below. 

A novelty of the present study is the calculation of the 
grain charge distribution, in parallel with the chemistry 
and the dynamics; the charge distribution is significant 
because ions recombine predominantly on the surfaces 
of negatively charged grains, at a rate which is determined 
by their opacity per H nucleus, ( 71 ^/ 7111 )^ 0 ^. The higher 
fractional abundance of free electrons at high gas densi¬ 
ties, owing to coagulation (cf. Fig.^, leads to an increase 
in the fraction of the grains which are negatively charged; 
see Fig. |5l The fraction of positively charged grains al¬ 
ways remains small. When there is no coagulation, the 
fraction of neutral grains remains above 10 % as the col¬ 
lapse proceeds, whereas, for the maximum rate of coagula¬ 
tion considered in Fig.01 approximately 95% of the grains 
are negatively charged by the time that nn = 10 ^ cm“^. 
Nonetheless, the opacity per H nucleus of the negatively 
charged grains, {ti^-/ 7111 ) 170 g, decreases as Vcrit increases, 
owing to the reduction in the total grain opacity per H 
nucleus as coagulation proceeds: see Fig. 0 l 3 . 

4.4. Nitrogen chemistry 

One problem which all models of the chemistry of prestel- 
lar cores have to solve is the variation of the fractional 
abundances of nitrogen-containing compounds with den¬ 
sity. Observations indicate that, whilst most molecules 
deplete on to grain surfaces at densities of a few times 
10^ cm“^, species containing nitrogen, such as NH 3 and 
N 2 H+, survive to densities which are an order of mag¬ 
nitude larger (e.g. Tafalla et al. 2004, Belloche & Andre 
2004, Bergin et al. 2002). Chemical modellers have at¬ 
tempted to explain this behaviour by assuming lower val¬ 
ues of the adsorption energy for N 2 than for CO and other 
abundant neutral molecules; this results in cosmic ray- 
induced or direct thermal desorption of N 2 being more ef- 
hcient than for CO. Aikawa et al. (1997, 2001) and Bergin 
et al. (2002) give = 750 K as the adsorption en¬ 
ergy of N 2 and = 960 K as the adsorption energy of 
CO, on CO ice. On the other hand, recent measurements 
(Oberg et al. 2004) indicate that the adsorption energies 
of N 2 to N 2 and CO to CO ice are = 790 K and 
Eqq = 855 K; we have adopted these values in our cal¬ 


culations but considered the implications of using larger 
values of (A^'q — E/f^^). We have concluded that the ap¬ 
proximately 10 % difference between the binding energies 
of CO and N 2 , as measured by Oberg et al. (2004), is insuf¬ 
ficient to account for the observed differences in fractional 
abundances, to which reference was made above. We note 
that, in the present context, it is the ratio between the 
rates of desorption of CO and N 2 which is significant and 
determined by the difference in their adsorption energies; 
the absolute rates are less crucial. 

An possible explanation of the observations is that 
there is a difference between the sticking coefficients of 
CO and one of the principal forms of gas-phase nitro¬ 
gen (N or N 2 ). Accordingly, we have carried out calcu¬ 
lations in which the sticking coefficient, S, of either N 
or N 2 was set equal to 0.1, whilst retaining S' = 1 for 
all other species; the results of the calculations in which 
S(N 2 ) = 0.1 are shown in Fig. 01 In Fig.0t; and Fig.01i, we 
compare the free-fall models, including coagulation and 
with S(N 2 ) = 1 or S(N 2 ) = 0.1, respectively. In the case 
of S(N 2 ) = 0 . 1 , the fractional abundance of N 2 H+, which 
traces N 2 , decreases rapidly at densities nn ^ 10 ® cm“®, 
which is roughly an order of magnitude higher than for 
the corresponding fall-off of HCO'*'. The fractional abun¬ 
dance of N 2 H+ attains values of several times 10 “^®, which 
are, in fact, larger than the observed values. However, the 
effect upon the evolution of the major ions, H+ and H;^, 
is not large, and N 2 H“'' never becomes a major ion. We 
conclude that a small sticking coefficient for N 2 could ex¬ 
plain the observations; but this assumption is clearly ad 
hoc and needs experimental verification. 

The other possibility is that the sticking coefficient of 
atomic nitrogen is small, S'(N) = 0.1. This assumption has 
the consequence that a substantial fraction of the nitrogen 
remains in the gas phase, but it does not result in a signif¬ 
icant enhancement of the abundance of N 2 H+, relative to 
the case in which all the sticking coefficients S' = 1. The 
reason for this result is to be found in the gas-phase chem¬ 
istry of nitrogen. The paths to N 2 H+ and NH 3 start with 
the neutral-neutral reactions N(OH, H)NO and NO(N, 
0 )N 2 (Pineau des Forets et al. 1990), whose rates depend 
on the fractions of N and OH in the gas phase. However, 
these reactions are so slow that an increase in the initial 
gas-phase abundance of N has little effect on the abun¬ 
dances of N 2 H''' or NH 3 on the timescale of the collapse. 
Only when the abundance of N 2 is increased directly, by 
reducing its sticking coefficient or its adsorption energy, 
are the abundances of N 2 H+ and NH 3 enhanced signifi¬ 
cantly. 

4.5. Validation of previous results 

In our previous calculations of complete depletion in pre- 
protostellar cores (Walmsley et al. 2004; Flower et al. 
2004), we computed abundances in steady-state, neglect¬ 
ing dynamical effects but including deuteration reactions. 
Although we do not currently consider deuteration reac- 
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tions in the phase of ‘free-fall’, it is possible to compare 
the fractional ionization of the gas and the abundances 
of the major ions, H"*" and over the range of den¬ 
sities considered previously (2 x 10 ® < uh < 2 x 10 ^ 
cm“®). In order to make the appropriate comparison 
with our previous work, we plot, in Fig. m n(H^)/nH 
from the present calculations, neglecting coagulation, and 
[ 71 ( 133 ) + 7i(H2D+) -I- n(HD 2 ) + 77(D[)’)]/nH from our pre¬ 
vious study (Walmsley et al. 2004). 

As may be seen from Fig. the agreement between 
the two calculations is good. Although the abundances of 
atomic ions such as Fe’*' and, especially, S+ are initially 
significant (i.e when tih ~ 10 ^ cm“®), recombination with 
electrons, on the surfaces of grains, neutralizes the atomic 
ions. By the time that a density of tih « 10® cm“® has 
been reached, the principal ions are H;)" and H’*', as pre¬ 
dicted by our earlier calculations, in steady-state. Thus, 
even when the collapse is as fast as in the case of ‘free-fall’, 
the assumption of a steady-state is a good approximation. 

5. Concluding remarks 

We have considered the early stages of isothermal col¬ 
lapse of a pre-protostellar core. The initial conditions were 
taken to be those appropriate to dense molecular clouds, 
with significant fractions of the heavy elements being de¬ 
pleted into the cores or the ice mantles of dust grains. The 
time for collapse is taken equal to the timescale which 
characterizes either free-fall or ambipolar diffusion; the 
latter is almost an order of magnitude larger than the for¬ 
mer. 

The processes of further freeze-out on to the interstel¬ 
lar grains and of coagulation of the grains were included 
when calculating the chemical evolution of the medium 
in the course of its collapse. Species containing heavy ele¬ 
ments freeze on to the grains in the early phases of the col¬ 
lapse. When ambipolar diffusion determines the timescale 
for collapse, freeze-out occurs at densities which are lower 
than the values which have been deduced from observa¬ 
tions of CO (cf. Bacmann et al. 2002, Tafalla et al. 2004); 
we conclude that the observational data favour a collapse 
timescale of order of the free-fall time, rather than the 
(larger) ambipolar diffusion time. In this case, coagula¬ 
tion is unlikely to be sufficiently rapid to influence the 
rate of freeze-out. 

The rate of grain coagulation is determined by the 
turbulent velocity. At;, of the medium. In the limit of 
Av « Cs, the rate of coagulation is negligible. On the 
other hand, when proceeding at close to its maximum rate, 
coagulation has important consequences for the degree of 
ionization and the ionic composition of the medium. At 
high gas densities, the major ions are H+ and H;^, and 
their relative abundance depends on the rate of coagu¬ 
lation, through its effect on the grain opacity per H nu¬ 
cleus, ( 7 ig/nH) 7 rag; this is because H+ ions recombine on 
the surfaces of negatively charged grains, whereas ions 
also recombine dissociatively with free electrons in the gas 
phase. If coagulation is slow or absent, n(H;^) > 77 (H+), 


whereas, if coagulation is rapid, 7 i(H[^) < ri(H+) at high 
densities. The observations of H 2 D+ and D 2 II+ in pre- 
protostellar objects suggest that Hg and its deuterated 
forms are likely to be the dominant ions, in which case 
coagulation is slow or absent. 

We consider that it is important to calculate the grain 
charge distribution, in parallel with the chemistry and the 
dynamics, because H"*' ions recombine predominantly on 
the surfaces of negatively charged grains. Grain coagu¬ 
lation results in an increase in the fraction of electrons 
which remain in the gas phase at high densities and to an 
increase in the fraction of the grains which are negatively 
charged; the fraction of positively charged grains is always 
small. 

Observations of pre-protostellar cores indicate that 
nitrogen-containing species, such as NH3 and N2H'*', sur¬ 
vive to densities which are an order of magnitude larger 
than those at which most molecules freeze on to grain sur¬ 
faces (Tafalla et al. 2004, Belloche & Andre 2004, Bergin et 
al. 2002). This differential depletion was believed to be due 
to the adsorption energy of N 2 being significantly lower 
than those of other molecules, such as CO (see Bergin & 
Langer 1997); but recent measurements of the adsorption 
energies of CO and N 2 (Oberg et al. 2004) do not support 
this interpretation. In order to solve this problem, a mech¬ 
anism must be found which results in nitrogen-containing 
species remaining in the gas phase at densities above 10 ® 
cm“®; the gas-phase chemistry must also be adequately 
described under these conditions. We have considered the 
consequences of reducing the sticking coefficient of N or 
N 2 to values of the order of 0.1. Only when S'(N 2 ) = 0.1, 
whilst retaining S' = 1 for all other species, are the abun¬ 
dances of both NH3 and N2H''' enhanced, relative to CO, 
during the course of the collapse. However, we recognize 
that the assumption of S(N 2 ) « 0.1 is ad hoc and needs 
experimental verification. 


Appendix A: The rate coefficient for grain 
coagulation 

We assume that the velocity distribution of the popula¬ 
tion of grains is Gaussian and consider a collision between 
two grains, 1 and 2. The x-component of the velocity dis¬ 
tribution of grain 1 is given by 


Af \ -hi 

71-2 


such that 


/ oo 

4>{Va;i)dVxi = 1 
-00 


(A.l) 


(A.2) 


Thus the probability that grain 1 has cartesian velocity 
components in the ranges —>■ + dT;^;^, Vy-^ Vy^ + 

dvy ^, ^ -I- dw^j is 


p{vxi,Vy^,v^^) dvxi dvy^ =(j){vxi) (l){'Vyi) 4>{uzi) dv^^ dVy^ dV:,^ 
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e dvxi dvy^ dv^^ (A.3) 

7r2 / 

where t;? = +v"^, +v‘^ and we have assumed the veloc- 

yi ^1 

ity distribution to be isotropic. An expression completely 
analogous to foil holds for grain 2. Following Margenau 
& Murphy (1956), it may be shown that the probability 
of grain 1 having an ^-component of velocity in the range 
Vxi —^ Vxi + dvxi and grain 2 having an ^-component of 
velocity in the range Vx^ —> Vx^ + dvx 2 is 


of the population of grains is charged, may be seen to be 
justified as follows. 

The grain population consists of neutral and singly 
charged grains; of the latter, negatively charged grains 
are by far the most abundant. Consider the long-range 
Coulomb interaction between two grains, each with a sin¬ 
gle negative charge; this repulsive potential is given by 

hCoul — ; 

r 


li tt2 2 

(j){vxC) di;,ci dvx^ = —e“ dvx 

7r2 

where 

1 _ 1 1 _ 2 
Ip ~ 


(A.4) 


where r is the separation of the grain centres; Vcoui has a 
maximum value of e^/(2ag), when the grains are in con¬ 
tact. The kinetic energy of a grain, in the centre of mass 
system, is TOgWg/4, where 


(A.5) 



when hi = /i 2 = h; Vx = Vx^ — Vx 2 is the relative x- 
component of velocity of the two grains. Analogous ex¬ 
pressions hold for the y- and z-components of velocity. 
Thus, the probability of finding the velocity components 
of the grains within the ranges specified above is 

p{uxi,Vy^,V2:J pivx2,Vy2,V2:^) dvx^ dvy^ dz;zi dvx2 dvy2 dv 
— e ” dvx dvy dvz 

7r2 / 

= P{vx,Vy,Vz) dvx dvy duz (A.6) 


is the mass of a grain (and mg/2 is the reduced mass, for 
collisions between identical grains); Vg is the relative col¬ 
lision speed. Taking Ug = 0.05/xm, which is the minimum 
(initial) value in the calculations reported in this paper, 
and pg = 2 g cm“^, we obtain nig = 1.0 x 10“^^ g. The rel- 
Z 2 ative kinetic energy of the grains is greater than or equal 
to the repulsive Coulomb potential when 

3 y e 

2T:pg) al 


where v = vi — V 2 is the relative velocity of the grains 1 
and 2. As the distributions are isotropic, 

dvx driy dvz = dr; (A.7) 

Defining 

V{v) = d'KP[Vx,Vy,Vz) 


i.e. when Vg ^ 10 cm s“^ for the parameters above. On 
the other hand, the Gaussian velocity distribution of the 
grains has a FWHM of Au = Cs/4 Ri 6000 cm s“^ in the 
calculations reported here; see Fig. ro Thus, only a very 
small section of the velocity distribution, with Vg ^ 10 
cm s“^ << At!, is affected significantly by the Coulomb 
interaction. The polarization potential between the grains 
is of shorter range and attractive; its asymptotic form is 


= 47r 



(A.8) 


Mpoi = - 


2r4 ’ 


we have the overall normalization condition 

P{v)v‘^ du = 1 (A.9) 

The rate coeSicient for coagulation in grain-grain col¬ 
lisions is 

r'^cxit 

(av) = dTTOg / vV{v)v'^dv (A.10) 

Jo 

where tr = diral is the cross section for collisions between 
identical grains; this is the expression for the geometrical 
cross section, assuming the grains to be spherical, which 
minimizes their surface area (and effective cross section) 
for a given volume. Clearly, non-spherical grains would 
have a larger cross section. The neglect of the Coulomb 
and polarization potentials, which arise when at least part 



where Cg is the polarizability of a grain of radius Ug. 
Taking into account this interaction would tend to en¬ 
hance the cross section for coagulation of a neutral and 
a charged grain; but arguments analogous to those imme¬ 
diately above show that the effect would be negligible for 
the velocity distribution that we have adopted. 

Following Chokshi et al. (1993), we assume that the 
grains stick (coagulate) if the relative collision speed 
V < fcrit- The critical speed, fcrit oc according to 

Chokshi et al. (1993); the constant of proportionality de¬ 
pends on the grain composition. 

Using (IA.8II . the integral in equation (ICTt can be 
evaluated, yielding 


{av) 


= a 


TT In 2 


2 


Av 



1 -h 2 In 2 



exp 
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= cr 0.9584 Ar; U - 


(l- 

1-b 1.3863 

exp 

-1.3863 

' f caind^ 

l 

\AvJ_ 



^ Ahhei dt 


where we have substituted = h^/2 and defined 
4 In 2 


= 


(Av) 


(A.12) 


where Av is the FWHM of the Gaussian distribu¬ 
tions 1131) and (Toil above, i.e. the width of the grain 
velocity distribution. 

Using equation mu, it may be shown that, when 

t^crit AlJ, 


iav) 


(2 In 2)3 


(Ar;)3 


are assumed to have a unique radius, Og, then 
, 3 


4agpg 


(B.3) 


where Mg = rjUYiPg is the mass of grains per unit volume of 
gas. Substituting for rj from equation (lOll . (Q becomes 


2 3 Tynn 

- 4 a 


Equating equations and we obtain 




- 0.5 




- a, 


- 0.5 


(B.4) 


(B.5) 


= 0.9209 

whereas, when Ucrit >> Av, 

(av) « 0.9584 aAv 

In practice, Ucrit and Av are comparable, and the exact 
expression mu was used to calculate the coagulation 
rate coefficient. 

In Fig. ro is plotted the variation of the rate coef¬ 
ficient for coagulation, (av), with the turbulent width. 
Ad, for Og = 0.05 /im. Also shown in the Figure are 
fractions of the adiabatic sound speed, Cs (evaluated at 
T = 10 K). The linear dependence of {av) on Ad, in the 
limit Ad << Dcut, is evident in this Figure. The rate co¬ 
efficient has a maximum near Ad = Cs/4, for the adopted 
value of Og and assuming that Dcrit(og) = 0.4 Og (Poppe 
& Blum 1997; broken curve in Fig. mi); turbulent widths 
much in excess of or much below this value result in neg¬ 
ligible coagulation. 


Appendix B: The initial value of Og 

Here we demonstrate that the value of Og required to 
reproduce a given value of rigTrOg depends only on the 
adopted size distribution and its limits, and not on the 
density of the grain material nor the dust:gas mass ra¬ 
tio. The grain size distribution is taken from Mathis et al. 
(1977): 

drig oc o“3'3(Jq^ 


Let om be the upper limit to the distribution and Um be 
the lower limit. Then, following Le Bourlot et al. (1995), 
the mean grain opacity is given by 




3 


- 0.5 


-, 0.5 


— a 


- 0.5 

m 


^0.5 


(B.l) 


where 


7 ] = 


lAmuG 

Pg 


(B.2) 


which is independent of pg and G. 

Acknowledgements. GdesF and DRF gratefully acknowledge 
support from the ‘Alliance’ programme, in 2004 and 2005. We 
are thankful to Ewine van Dishoeck and Helen Fraser for com¬ 
municating results of Oberg et al. (2004), in advance of their 
publication, and to an anonymous referee for a thoughtful and 
detailed report. 

References 

Aikawa, Y., Miyama, S. M., Nakano, T., Umebayashi, T. 1996, 
ApJ, 467, 684 

Aikawa, Y., Umebayashi, T., Nakano, T., Miyama, S. M. 1997, 
ApJ, 486, L51 

Aikawa, Y., Ohashi, N., Inutsuka, S., Herbst, E., Takakuwa, S. 
2001, ApJ, 552, 639 

Aikawa, Y., Ohashi, N., Herbst, E. 2003, ApJ, 593, 906 
Aikawa, Y., Herbst, E., Roberts, H., Caselli, P. 2004, 
astro-ph 0410582 

Anders, E., Grevesse, N. 1989, Geochim. Gosmochim. Acta, 53, 
197 

Bacmann, A., Lefloch, B., Ceccarelli, C., Castets, A., 
Steinacker, J., Loinard, L. 2002, A&A, 389, L6 
Belloche, A., Andre, P. 2004, A&A, 419, 35 
Bergin, E. A., Langer, W. D. 1997, ApJ, 486, 316 
Bergin, E.A., Alves, J., Huard, T., Lada, C.J. 2002, ApJ, 570, 
LlOl 

Bianchi, S., Goncalves, J., Albrecht, M. et al. 2003, A&A, 399, 
L43 

Caselli, P., van der Tak, F. F. S., Ceccarelli, C., Bacmann, A. 
2003, A&A, 403, L37 

Chokshi, A., Tielens, A. G. G. M., Hollenbach, D. 1993, ApJ, 
407, 806 

Ciolek, G. E., Basu, S. 2000, ApJ, 529, 925 
Draine, B. T., Sutin, B. 1987, ApJ, 320, 803 
Flower, D. R., Pineau des Forets, G. 2003a, MNRAS, 341, 1272 
Flower, D. R., Pineau des Forets, G. 2003b, MNRAS, 343, 390 
Flower, D. R., Pineau des Forets, G., Walmsley, C. M. 2004, 
A&A, 427, 887 

Geppert, W.D., Thomas, R., Semaniak, J., et al. 2004, ApJ 
609, 459 

Gibb, E. L., Whittet, D. C. B., Schutte, W. A., Boogert, A. C. 
A., Chiar, J. E., Ehrenfreund, P., Gerakines, P. A., Keane, 
J. V., Tielens, A. G. G. M., van Dishoeck, E. F., Kerkhof, 
O. 2000, ApJ, 536, 347 



















D.R. Flower et al.: Freeze-out and coagulation in pre-protostellar collapse 


11 


Hasegawa, T. L, Herbst, E. 1993, MNRAS, 261, 83 
Kramer, C., Richer, J., Mookerjea, B., Alves, J., Lada C. 2003, 
A&A, 399, 1073 

Le Bourlot, J., Pineau des Forets, G., Roueff, E., Flower, D. 
R. 1995, A&A, 302, 870 

Lee, J.-E., Bergin, E.A., Evans, N.J. 2004, ApJ, 617, 360 
Leger, A., Jura, M., Omont, A. 1985, A&A, 144, 147 
Margenau, H., Murphy, G. M. 1956, The Mathematics of 
Physics and Ghemistry (Van Nostrand: Princeton, N.J.) 
Mathis, J. S., Rumpl, W., Nordsieck, K. H. 1977, ApJ, 217, 
425 

Nakano, T., Umebayashi, T. 1980, PASJ, 32, 613 

Oberg, K. L, van Broekhuizen, F., Fraser, H. J., Bisschop, S. 

E. van Dishoeck, E. F., Schlemmer, S. 2005 ApJ, 621, L33 
Ossenkopf, V. 1993, A&A, 280, 617 
Ossenkopf, V., Henning, Th. 1994, A&A, 291, 943 
Osterbrock, D. E. 1961, ApJ, 134, 270 

Pineau des Forets, G., Roueff, E., Flower, D. R. 1990, MNRAS, 
244, 668 

Poppe, T., Blum, J. 1997, Adv. Sp. Res., 20, 1595 
Roberts, H., Herbst, E., Millar, T.J. 2003, ApJ, 591, L41 
Savage, B. D., Sembach, K. R. 1996, ARA&A, 34, 279 
Shu, F. H., Adams, F. G., Lizano, S. 1987, ARA&A, 25, 23 
Sofia, U. J., Meyer, D. M. 2001, ApJ, 554, L221 
Spitzer, L. 1978, Physical Processes in the Interstellar Medium 
(Wiley: New York) 

Suttner, G., Yorke, H. W. 2001, ApJ, 551, 461 
Tafalla, M., Myers, P. G., Caselli, P., Walmsley, C. M. 2004, 
A&A, 416, 191 

Vastel, C., Phillips, T. G., Yoshida, H. 2004, ApJ, 606, L127 
Wakelam, V., Gaselli, P., Geccarelli, G., Herbst, E., Castets, 
A. 2004, A&A, 422, 159 

Walmsley, C. M., Flower, D. R., Pineau des Forets, G. 2004, 
A&A, 418, 1035 





Fig. 1. The fractional abundances of (a) selected neutrals 
and (b) the major ions, as functions of the evolution time, 
t. Results are shown for the limiting case of no coagulation 
and assuming that the timescale for collapse is determined 
by the free-fall time, rg. The corresponding density pro¬ 
file, nH(t), is plotted in panel (c). 
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Fig. 3. The variation of the grain radius, Og, as a function 
of the gas density. Results are shown for the limiting case 
of no coagulation (ucrit = 0) and for Ucrit(ag) = 0.1 Og 
(Chokshi et al. 1993) and Ucrit(ag) = 0.4 (Poppe & 
Blum 1997); Av = Cs/4, where Av is the turbulent width 
and Cs is the adiabatic sound speed. The timescale for 
collapse is determined by the free-fall time, rg. 


Fig. 2. The fractional abundances of the exemplary neu¬ 
trals, N 2 and CO, as functions of the gas density, nn- 
Results are shown (a) for the limiting case of no coagu¬ 
lation, with the timescale for collapse determined by ei¬ 
ther the free-fall time, rg, or the ambipolar diffusion time, 
Tad; (b) including coagulation, with Ucrit(ag) = 0.1 Ug 
(Chokshi et al. 1993; full curves) or Ucrit(ag) = 0.4 Og^^® 
(Poppe & Blum 1997; broken curves). Av = Cs/4, where 
Av is the turbulent width and Cg is the adiabatic sound 
speed. The timescale for collapse is determined by the 
free-fall time. 
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Fig. 4. The fractional abundances of the specified gas- 
phase ions and the electrons, as functions of the gas den¬ 
sity. Results are shown for (a) the limiting case of no coag¬ 
ulation (ucrit = 0), (b) Ucrit(ag) =0.1 Og [equation Q; 
Chokshi et al. (1993)], and (c) Ucrit(ag) = 0.4 Og (Poppe 
& Blum 1997); Av = Cs/4, where Av is the turbulent 
width and Cg is the adiabatic sound speed. The results 
in panel (d) were obtained on reducing the value of the 
sticking coefficient for molecular nitrogen to S'(N 2 ) = 0.1 
and should be compared with those in panel (c), where 
S'(N 2 ) = 1. In all cases, the timescale for collapse is deter¬ 
mined by the free-fall time, rg. 


Fig. 5. (a) The fraction and (b) the opacity per H nucleus, 
{n~/nYi)Trag, of negatively charged grains, as functions of 
the gas density, nn. Results are shown for the limiting 
case of no coagulation (ucrit = 0), for i;crit(ag) = 0.1 Og 
[equation Chokshi et al. (1993)], and for i;crit(ag) = 
0.4 (Poppe & Blum 1997); Av = Cs/4, where Av is 

the turbulent width and Cg is the adiabatic sound speed. 
The timescale for collapse is determined by the free-fall 
time, Tff. 
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Fig. A.l. The variation of the rate coefficient for coagu¬ 
lation, (crn), with the turbulent width, Av, for Og = 0.05 
/rm and Ucrit(ag) = 0.1 (Chokshi et al. 1993; full 

curve) or t'crit(ag) = 0.4 Og (Poppe & Blum 1997; bro¬ 
ken curve). Also shown are fractions of the adiabatic sound 
speed, Cs, evaluated at T = 10 K. 


Fig. 6 . A comparison of the fractional abundances of the 
principal ions, (a) as computed by Walmsley et al. (2004), 
assuming steady-state, and (b) as predicted by the present 
dynamical model; Og = 0.05 /rm and, in this case only, = 
3 X 10“^^ s“^, in order to be able to compare directly with 
Walmsley et al. (2004). [“Hg ” denotes ^(Hg )-|-n(H 2 D+)-|- 
n(HD+)-fn(D+).] 





